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ABSTRACT 

We obtain a solution for the hydrodynamic outflow of the polytropic gas 
from the gravitating center, in presence of the uniform Dark Energy (DE). The 
antigravity of DE is enlightening the outflow and make the outflow possible at 
smaller initial temperature, at the same density. The main property of the wind 
in presence of DE is its unlimited acceleration after passing the critical point. 
In application of this solution to the winds from galaxy clusters we suggest that 
collision of the strongly accelerated wind with another galaxy cluster, or with 
another galactic cluster wind could lead to the formation of a highest energy 
cosmic rays. 

Subject headings: dark energy — - galaxy clusters — -galactic wind 

1. Introduction 

It was shown by Chernin (2001, 2008) that outer parts of galaxy clusters (GC) may 
be under strong influence of the dark energy (DE), discovered by observations of SN la at 
redshift z < 1 (Riess et al., 1998; Perlmutter et al., 1999), and in the spectrum of fluctua- 
tions of the cosmic microwave background radiation (CMB), see e.g. Spergel et al. (2003), 
Tegmark et al. (2004). Equilibrium solutions for polytropic configurations in presence of 
DE have been obtained in papers of Balaguera-Antolinez et al. (2006, 2007), and Merafina 
et al. (2012). The hot gas in the galactic clusters may flow outside due to high thermal 
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pressure, and in the outer parts of the cluster the presence of a dark energy (DE) facilitates 
the outflow. 

Here we obtain a solution of hydrodynamic equations for the winds from galactic clus- 
ters in presence of DE. We generalize the solution for the outflows from the gravitating 
body, obtained for solar and stellar winds by Stanyukovich (1955) and Parker (1963), to the 
presence of DE. It implies significant changes in the structure of solutions describing galactic 
winds. 



2. Newtonian approximation in description galactic winds in presence of DE 

A transition to the Newtonian limit, where DE is described by the antigravity force in 
vacuum was done by Chernin (2008). In the Newtonian approximation, in presence of DE, 
we have the following hydrodynamic Euler equation for the spherically symmetric outflow 
in the gravitational field of matter and DE 

dv dP ( Gm m Ac 2 r\ ( Gm m 8nGp\r\ 
pv— + — = -p[ — — = -p I — ] . 



dr dr \ r 2 3 / \ r 2 3 / 

Here p and P are a matter density and pressure, respectively, m m is the mass of the matter 
inside the radius r. We use here DE in the form of the Einstein cosmological constant 
A. Newtonian gravitational potentials produced by matter $ 9 , and $a by DE, satisfy the 
Poisson equations 



A$ A = -87rGp A , A$ g = 4nGp, p A = —. (2) 



Ac 2 

We consider, for simplicity, the outflow in the field of a constant mass (like in stellar wind) 
m m = M. The Eq. (TjQ) in this case is written as 



dv dP (GM Ac 2 r\ ( GM %TtGp k r\ 

The Eq. (pQ) should be solved together with the continuity equation in the form 



A-Kpvr 2 = M, (4) 

where M is the constant mass flux from the cluster. We consider polytropic equation of 
state, where pressure P, and sound speed c s are defined as 
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Introduce nondimensional variables as 



v _ c, _ r GM 



C. 



f* = c*, (6) 



1 7 



i 



In non-dimensional variables the equation (j3J) is written as 

_dt> 2 _ dc s 1 , n . Ac 2 r 2 

^T" + + ^ - Ar = 0, A = — ^. 7 

The continuity equation in non-dimensional form is written as 

2 /t> 

pt>r = 771, Cs~ fr = 777., 777 = — M* = 47rp*'U*7\ ) .. (8) 
It follows from that 

(i5 2 dc, tip dv df 

4- = ~, 4- + — + 2— = 0. 9 

p 7 — 1 c s p i> r 

Using ([9]) we may write the equation of motion (JTj) in the form 



cfti v 2c 2 - I + Af 2 



(10) 



or r t> 2 — c 2 

The only physically relevant solutions are those which pass smoothly the sonic point v — c s , 
being a singular point of the Eq. (10), with 

v = 5 s , 2d 2 -- + \r 2 = (11) 
r 

where r = r c , v = v c , c s = c sc . Choosing c* = c sc , we obtain in the critical point 

v c = c sc = 1, 2 - 1 + Af 2 = 0. (12) 
With this choice of the scaling paraneters, we have from (jSj) 
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m = f 2 . (13) 

The relation f|T2|) determining the dependence r c (A) in the solution for the galactic wind and 
accretion, in presence of DE, is presented in Fig.l. The Eq.flT]) for the polytropic flow has a 
Bernoulli integral as 




Fig. 1. — The functions r c (A) (full curve), according to ( |T2l) . and f max (X) (dashed curve) 
according to ( l22p . It is clear, that the critical radius of the flow r c is always inside the radius 
of the extremum of the total gravitational potential r max . 



v 2 c 2 1 Ar 2 , _ 2 /m V ' / r 2 x ' 



H ^-r - ~ 7T = hi c = — = U% ■ ( 14 ) 

2 7 — 1 r 2 \vr z / \vr z J 

The dimensional Bernoulli integral H = hc 2 c . The Bernoulli integral is determined through 
the parameters of the critical point, with account of f lT2"j) . as 

2(7-1) f c 2 2(7-1) 2\f c J 1 ' 

The dependence h(X) for different polytropic powers 7 is given in Fig. 2. 
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Fig. 2. — The function h(X) for 7 = | (full curve); 7 = | (dashed curve); 7 
curve), according to relations f lT2|) .f lT5|) . 
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The stationary solution for the wind is determined by two integrals: constant mass flux 
M, and energy (Bernoulli) integral H. In absence of DE we obtain the known relations 



1 , 5-3 7 
-> h — — 

2' 2(7-1) 



(16) 



At small A we have from (IT2]),([T5 



A 



5-37 3 
— ~A 



r c = 0.5 , h 

16' 2(7-1) 



(17) 



At large A — >■ 00 it follows from (fT2|) r c — > r, 
around r coo in the form 



A 1 / 3 . Making expansion in (|T2|) 



I = X 1/3 + e, 



we obtain from (fT2|) . f lT5|) 
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£ = -, r c = — ^_ fe = 5 37 - -^ + 2= 7— — r — -A 1/3 at A ->■ oo. (18) 
3' ^ 1/3 + | 2(7-1) 2 2(7-1) 2 v ; 

In the outflow from the physically relevant quasi-stationary object the antigravity from DE 
should be less the the gravitational force on the outer boundary, which we define at r = r*. 
Therefore the value of A is restricted by the relation (see e.g. Bisnovatyi-Kogan and Chernin, 
2012) 

Ac 2 AnM 

2 ^ = ^g < ^^t (19) 

In non-dimensional variables this restriction, with account of ©,© is written as 

16tt 2 

A < — = 17.55 = \ lim . (20) 



It is reasonable to consider only the values of A smaller than Xu m . It follows from (fl2]) . 

that f c is monotonically decreasing with increasing A. For A = Xu m = 17.55 we obtain 
?c = ^c,Um ~ 0.29. The effective gravitational potential $ is formed by the gravity of the 
central body, and antigravity of DE 

1 Xr 2 

$ = . 21 

f 2 v ; 

To overcome the gravity of the central body, the value of h should exceed the maximum 
value of the gravitational potential, defined by the extremum of $ 

h > $ (f ) = --A 1/3 f = A~ 1/3 (22) 

11 — ^maxy'max) ^ > 'max /y ■ \^^J 

The function f max (X) is represented in Fig.l, it is always f max > f c . So, in presence of DE 
the outflow of the gas from the cluster to the infinity is possible even at the negative values 
of h. In absence of DE the non-negative value of h, and the outflow are possible only at 

7<§. 



3. Solutions of the galactic wind equation in presence of DE 
3.1. Analytical solutions 

Analytical solutions exist only in absence of DE, A = 0. In this case we have from 

i a-dna 
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f-i fe-* I ^ *- 5 ~ 3T (23) 
rc "2' / '" 2 + 7 -l f~ 2(7-1)' ( j 

For 7 = | the constant h = for the critical solution. At /i = and 7 = | there is a whole 

family of solutions with arbitrary constant Mach number Ma=^- in all space. It is more 

convenient to write this solution in dimensional variables, with Bernoulli integral 

TT v 2 3c 2 GM . . 

H = — + — £ = 0. 24 

2 2 r K J 

There is an exact solution in the form 

2 _ 2GM Ma 2 2 _ 2GM 1 

U ~ r 3 + Ma 2 ' ° s ~ r 3 + Ma 2 ' ( ^ 

At Ma=l this solution corresponds to the critical solutions for 7 = §• 

Another analytic solution takes place at 7 = 1.5, with the non-dimensional Bernoulli 
integral h = 0.5. Using from (j!4p . for 7 = 1.5, r c = 0.5, 



^2 \ 7-1 



7,2 I r c 



Avr 2 ' 



we obtain 

v 2 1 1 1 

/l — ~ — h 



2 rVti r 2 
This nonlinear equation has two solutions 



1.0 = 1; 2. f=-= = . (26) 

+ V^)(l + i3) 

There, the first solution corresponds to the wind, and the second one corresponds to the 
accretion case, with negative velocity; v in ff26l) represents the absolute value of this variable. 



3.2. Numerical solution 

To obtain a physically relevant critical solution of (iID]L with c 2 from (1141) . we obtain 
expansion in the critical point with v 2 = c 2 = 1, in the form 



v — 1 + a(r — r c ), ai 
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2 7 - 1 12 I 1 Af c 2 / 1 Xr 2 c 

a = ! a/2 H 1 --7(2 

r c 7 + 1 f c 7 + 1 V 4f C 2 V 4f: c 2 

Here «i corresponds to the wind solution, and a 2 is related to the case of accretion where v 
define the absolute value. At A = we have a well known expansion with 



a i 



Oil 



7 + 1 



7 + 1 
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(7-1) 
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It follows from the expansion ( 1271) . that physically relevant solutions exist only with positive 
value under the square root. It give the restriction for the value of 7 as a function of A in 
the form 



9 1 J_ , Af| 
< Z 4f e 2 
7 ^ 7maz - l ^2 • 

Z 4f c 2 

At A = 32, f c = 0.25 the limiting value 7 maa; goes to 00, so that at A > 32 the wind solutions 
exist formally for all polytropic powers 7. The dependence 7 ma x(A) is given in Fig.3. 

The critical solutions of the equation (flQj) . with account of (|T4l . are presented in Figs. 4-6 
for different values of 7 and A. Both wind and accretion solutions are presented. 



4. Discussion 



It is clear that the presence of DE tends to help the outflow of the hot gas from the 
gravitating object, as well as to the escape of rapidly moving galaxies (Chernin et al, 2013). 
Here we have obtained the solution for outflow in presence of DE, which generalize the well- 
known solution for the polytropic solar (stellar) wind. Presently the DE density exceed the 
density of the dark matter, and, even more, the density of the barionic matter. The clusters 
which outer radius is approaching the zero gravity radius, may not only loose galaxies, which 
join the process of Hubble expansion, but also may loose the hot gas from the outer parts 
of the cluster. Let us consider outer parts of the Coma cluster at radius Rc = 15 Mpc, with 
the mass inside Mq = 5 • 1O 15 M , from Chernin et al. (2013). For the present value of 
Pa = 0.71 • 10~ 29 g/cm 3 , supposing that Rc = t* is the critical radius of the wind, we obtain 
from (2), (7), the nondimensional constant A as 
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Fig. 3. — The functions 7 ma: r(A). The smooth solution for the outflow from the gravitating 
center, in presence of DE, is possible only at 7 < 7 ma:r . 



It corresponds to the temperature about T 6 ■ 10 7 K, A;T 5 keV. Observations of the 
hot gas distribution in the Coma cluster (Watanabe et al., 1999) on ASCA satellite have 
shown a presence of hot region with kT — 1 1 — 14 keV, and more extended cool region with 
kT — 5 ± 1 keV, what is in good accordance with our choice of parameters. 

Wind solutions for A=0; 0.58; 1.1 are presented in Fig. 7. The solution with A=0.58 is 
the closest to the description of the outflow from Coma cluster. The density of the gas in the 
vicinity of r = r c is very small, so the flow may be considered as adiabatic (polytropic) with 
the power 7=5/3. Without DE such gas flow is inefficient, its velocity is decreasing ~ 1/y/r, 
see Eq. (25). In presence of DE the wind velocity is increasing 2 times at the distance of 
~ 5r c ~ 75 Mpc from Coma. 

After quitting the cluster the gas is moving with acceleration, acting as a snowplough for 
the intergalactic gas. The shell of matter, forming in such a way, may reach a high velocity, 



Fig. 4. — The integral curves of the equations (flOj) . (JHJ), for 7 = 4/3 and A = 0, r c =0.5 
(dashed curves); A = 1.10, r c = 0.45 (dash-dot-dot curves); and A = 5.13, r c = 0.37 (full 
curves). Wind solutions correspond to curves with increasing velocity at large radius. The 
curves with decreasing velocities correspond to the accretion solution with negative v , so 
that its absolute value is presented. 

exceeding considerably the speed of galaxies in cluster. If the shell meets another cluster, or 
another shell moving towards, the collision of such flows may induce a particle acceleration. 
Due to high speed, large sizes, and low density such collisions may create cosmic rays of 
the highest possible energy (EHECR). We may expect the largest effect when two clusters 
move to each other. The influence of DE is decreasing with with a red shift, therefore the 
acceleration of EHECR in this model should take place in the periphery, or between, the 
closest rich galaxy clusters. 
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Fig. 7. — The integral curves of the equations ffTUj) . ffT4|) for the wind solution, at 7 = 5/3 
and A = 0, r c =0.5 (dashed curve); A = 1.10, r c = 0.45 (dash-dot-dot curves); and A = 0.58, 
r c = 0.47 (full curves). 



